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Abstract 



The prompt production of charmonium Xc and J/ip states is studied in proton-proton 
collisions at a centre-of-mass energy of y/s = 7 TeV at the Large Hadron Collider. The 
Xc and J/ip mesons are identified through their decays Xc —> J/^l and J/ip — > fi + fj,~ 
using 36 pb _1 of data collected by the LHCb detector in 2010. The ratio of the prompt 
production cross-sections for Xc and J/ip, cr{Xc J/i J l)/< J (J/i J ), is determined as a 
function of the J/ip transverse momentum in the range 2 < p'^ < 15 GeV/c. The 
results are in excellent agreement with next-to-leading order non-relativistic expectations 
and show a significant discrepancy compared with the colour singlet model prediction at 
leading order, especially in the low region. 
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1. Introduction 

The study of charmonium production provides an important test of the underlying 
mechanisms described by Quantum Chromodynamics (QCD). At the centre-of-mass en- 
ergies of proton-proton collisions at the Large Hadron Collider, cc pairs are expected to 
be produced predominantly via Leading Order (LO) gluon-gluon interactions, followed by 
the formation of bound charmonium states. The former can be calculated using pertur- 
bative QCD and the latter is described by non-perturbative models. Other, more recent, 
approaches make use of non-relativistic QCD factorization (NRQCD), which assumes the 
cc pair to be a combination of colour-singlet and colour-octet states as it evolves towards 
the final bound system via the exchange of soft gluons [lj . The fraction of J/i/j produced 
through the radiative decay of Xc states is an important test of both the colour-singlet and 
colour-octet production mechanisms. In addition, knowledge of this fraction is required 
for the measurement of the J/0 polarisation, since the predicted polarisation is different 
for J/ip mesons coming from the radiative decay of Xc state compared to those that are 
directly produced. 

In this paper, we report the measurement of the ratio of the cross-sections for the 
production of P-wave charmonia % c j(lP), with J = 0, 1, 2, to the production of J/ip in 
promptly produced charmonium. The ratio is measured as a function of the J/ip transverse 
momentum in the range 2<p^ < 15 GeV/c and in the rapidity range 2.0 <y J ^ <4.5. 
Throughout the paper we refer to the collection of Xej(l-P) states as Xc- The Xc and J/i/j 
candidates are reconstructed through their respective decays Xc — > J/07 an d J/0 — * 
using a data sample corresponding to an integrated luminosity of 36 pb _1 collected during 
2010. Prompt (non-prompt) production refers to charmonium states produced at the 
interaction point (in the decay of 6-hadrons); direct production refers to prompt J/0 
mesons that are not decay products of an intermediate resonant state, such as the ip(2S). 
The measurements are complementary to the measurements of the J/0 production cross- 
section [2] and the ratio of the prompt Xc production cross-sections for the J = 1 and J = 2 
spin states [3] , and extend the coverage with respect to previous experiments jU [5] . 

2. LHCb detector and selection requirements 

The LHCb detector [6] is a single-arm forward spectrometer with a pseudo-rapidity 
range 2<r]<5. The detector consists of a silicon vertex detector, a dipole magnet, a 
tracking system, two ring-imaging Cherenkov (RICH) detectors, a calorimeter system 
and a muon system. 

Of particular importance in this measurement are the calorimeter and muon systems. 
The calorimeter system consists of a scintillating pad detector (SPD) and a pre-shower 
system, followed by electromagnetic (ECAL) and hadron calorimeters. The SPD and pre- 
shower are designed to distinguish between signals from photons and electrons. The ECAL 
is constructed from scintillating tiles interleaved with lead tiles. Muons are identified using 
hits in muon chambers interleaved with iron filters. 

The signal simulation sample used for this analysis was generated using the Pythia 6.4 
generator [7] configured with the parameters detailed in Ref. [8]. The EvtGen [9], 
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Photos [TU] and Geant4 [TT] packages were used to decay unstable particles, generate 
QED radiative corrections and simulate interactions in the detector, respectively. The 
sample consists of events in which at least one J/ip — » / u + /i~ decay takes place with no 
constraint on the production mechanism. 

The trigger consists of a hardware stage followed by a software stage, which applies a 
full event reconstruction. For this analysis, events are selected which have been triggered 
by a pair of oppositely charged muon candidates, where either one of the muons has 
a transverse momentum > 1-8 GeV/c or one of the pair has pt > 0.56 GeV/c and the 
other has pt > 0.48 GeV/c. The invariant mass of the candidates is required to be greater 
than 2.9 GeV/c 2 . The photons are not involved in the trigger decision for this analysis. 

Photons are reconstructed using the electromagnetic calorimeter and identified using a 
likelihood-based estimator, CL 7 , constructed from variables that rely on calorimeter and 
tracking information. For example, in order to reduce the electron background, candidate 
photon clusters are required not to be matched to the trajectory of a track extrapolated 
from the tracking system to the cluster position in the calorimeter. For each photon 
candidate a value of CL 7 , with a range between (background-like) and 1 (signal- like), 
is calculated based on simulated signal and background samples. 

The photons are classified as one of two types: those that have converted to electrons 
in the material after the dipole magnet and those that have not. Converted photons 
are identified as clusters in the ECAL with correlated activity in the SPD. In order to 
account for the different energy resolutions of the two types of photons, the analysis 
is performed separately for converted and non-converted photons and the results are 
combined. Photons that convert before the magnet require a different analysis strategy 
and are not considered here. The photons used to reconstruct the Xc candidates are 
required to have a transverse momentum p T > 650 MeV/ c, a momentum p 1 > 5 GeV/ c 
and CL 7 > 0.5; the efficiency of the CL 7 cut for photons from Xc decays is 72%. 

All J/ip candidates are reconstructed using the decay J/ip — > n + n~ . The muon and 
J/ip identification criteria are identical to those used in Ref. [2J: each track must be 
identified as a muon with px > 700MeV/c and have a track fit % 2 /ndf <4, where ndf is 
the number of degrees of freedom. The two muons must originate from a vertex with a 
probability of the vertex fit greater than 0.005. In addition, the yU + /i~ invariant mass is 
required to be in the range 3062 — 3120 MeV/c 2 . The \ c candidates are formed from the 
selected J/ip candidates and photons. 

The non-prompt J/ip contribution arising from 6-hadron decays is taken from Ref. [2] . 
For the Xc candidates, the J/ip pseudo-decay time, t z , is used to reduce the contribution 
from non-prompt decays, by requiring t z = {zj/^ — zpy)Mji^ / p z < 0.1 ps, where Mj/^ is 
the reconstructed dimuon invariant mass, zm — zpv is the z separation of the reconstructed 
production (primary) and decay vertices of the dimuon, and p z is the z-component of the 
dimuon momentum. The z-axis is parallel to the beam line in the centre-of-mass frame. 
Simulation studies show that, with this requirement applied, the remaining fraction of Xc 
from 6-hadron decays is about 0.1%. This introduces an uncertainty much smaller than 
any of the other systematic or statistical uncertainties evaluated in this analysis and is 
not considered further. 



2 



The distributions of the n + fi~ mass of selected J/0 candidates and the mass difference, 
AM = M (/i + yU~ 7) — M (/i + yU~), of the selected Xc candidates for the converted and non- 
converted samples are shown in Fig. [I] The total number of prompt J/0 candidates 
observed in the data is ~ 2.6 million. The fit procedure to extract the three Xc signal yields 
using Gaussian functions and one common function for the combinatorial background is 
discussed in Ref. [3]. The total number of Xco, Xd an d Xc2 candidates observed are 823, 
38 630 and 26114 respectively. Since the Xco J/07 branching fraction is ~ 30 (17) 
times smaller than that of the xa (Xc2), the yield of Xco is small as expected [12]. 

3. Determination of the cross-section ratio 

The main contributions to the production of prompt J/0 arise from direct production 
and from the feed-down processes Xc — > J/07 and ip{2S) — > J/ipX where X refers to any 
final state. The cross-section ratio for the production of prompt J/0 from Xc — > J/07 
decays compared to all prompt J/0 can be expressed in terms of the three Xcj ( J = 0, 1, 2) 
signal yields, N XcJ , and the prompt J/0 yield, Nj/^, as 

°(x c -+ J/07) q(xc -» J/07) 

#) ~ a^(J/ij) + (7(0(25-) -> J/0X) + a( X c J/07) 



with 



and 



Njfy R 



2S 



J=2 

E 

J=0 
J=2 

+£ 

j=o 



N XcJ 



-sel 



A 7 , 



XcJ 



~XcJ ,Xc7 



-7 



-sel 



c dir 

,XcJ 



,dir 

,XcJ 



i? 



25 



(1) 



R 



25 



1 + /: 



25 



1 + /: 



,25 
2 ^ ,dir 



(2) 



./': 



25 



<t(i/j{2S) J/tpX) 



a 



dir 



WO 



(3) 



The total prompt Xc — > J/^7 cross-section is a(xc — > J/07) = Sj=o a Xcj'^(XcJ J/07) 
where a XcJ is the production cross-section for each Xcj state and B(xcj J/07) i s 
the corresponding branching fraction. The cross-section ratio /2s is used to link the 
prompt ip(2S) contribution to the direct J/0 contribution and -R25 takes into account 
their efficiencies. The combination of the trigger, reconstruction and selection efficiencies 
for direct J/0, for J/0 from 0(25") decay, and for J/0 from Xc — > J/07 decay are e^L, 
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M(h + h y) - M(nV) [MeV/c 2 ] 
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M(n + H7) - M(^V) [MeV/c 2 ] 

Figure 1: (a) Invariant mass of the /i + ' [T~ pair for selected J/ip candidates. The solid red curve 
corresponds to the signal and the background is shown as a dashed purple curve, (b) and 
(c) show the AM = M (fi + fj,~ j) — M(fi + /i~) distributions of selected Xc candidates with (b) 
converted and (c) non-converted photons. The upper solid blue curve corresponds to the overall 
fit function described in Ref. [3]. The lower solid curves correspond to the fitted XcO, Xd and \c2 
contributions from left to right, respectively (the XcO peak is barely visible). The background 
distribution is shown as a dashed purple curve. 
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e^, and e*/^ respectively. The efficiency to reconstruct and select a photon from a 
Xc — > J/07 decay, once the J/0 is already selected, is e* cj and the efficiency for the 
subsequent selection of the Xcj is e* c j J . 

The efficiency terms in Eq. (1) are determined using simulated events and are partly 
validated with control channels in the data. The results for the efficiency ratios e 2 ^ / 6jJ^ , 



,XcJ 



j/ii / e j/ii> 



and the product e^ cj e^ are discussed in Sect. |I 
The prompt Nj/^ and N XcJ yields are determined in bins of in the range 

2<pf'<15 GeV/c using the methods described in Refs. [2] and [3] respectively. In 
Ref. [2] a smaller data sample is used to determine the non-prompt J/0 fractions in 
bins of and rapidity. These results are applied to the present J/0 sample without 
repeating the full analysis. 



4. Efficiencies 

The efficiencies to reconstruct and select J/0 and Xc candidates are taken from simu- 
lation. The efficiency ratio £j/tyi/ e j/y> is consistent with unity for all p^ bins; hence, R 2 s is 
set equal to 1 in Eq. (|2). The ratio of efficiencies e^/e^ and the product of efficiencies 
e * cje sei f° r the Xci and Xc2 states are shown in Fig. [2j In general these efficiencies are 
the same for the two states, except at low where the reconstruction and detection 
efficiencies for Xc2 are significantly larger than for Xci- This difference arises from the 
effect of the requirement p^ > 650 MeV/ c which results in more photons surviving from 
Xc2 decays than from Xci decays. 
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Figure 2: (a) Ratio of the reconstruction and selection efficiency for direct J/0 compared to J/0 
from Xc decays, e^/e^, and (b) the photon reconstruction and selection efficiency multiplied 
by the Xc selection efficiency, c* cJ e^i , obtained from simulation. The efficiencies are presented 
separately for the Xci (red triangles) and Xc2 (inverted blue triangles) states, and as a function 
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The photon detection efficiency obtained using simulation is validated using candidate 
B + — > J/tpK + and B + — > XcK + (including charge conjugate) decays selected from the 
same data set as the prompt J/ip and Xc candidates. The efficiency to reconstruct and 
select a photon from a Xc in B + — > XcK + decays, e 7 , is evaluated using 



where N B +^ XcK + and N B +^ J / 1 p K + are the measured yields of B + — > XcK + and B + — > 
J/ipK + and B are the known branching fractions. The factor R e = 1.04 ±0.02 is obtained 
from simulation and takes into account any differences in the acceptance, trigger, selection 
and reconstruction efficiencies of the K , J/ip , Xc (except the photon detection efficiency) 
and B + in B + — > J/tpK + and B + — > XcK + decays. All branching fractions are taken from 
Ref. [12]. The B + J/ipK + branching fraction is B(B + J/^K + ) = (1.013 ± 0.034) x 
10~ 3 . The dominant process for B + — >• XcK + — )■ J/ip^K + decays is via the xa state, 
with branching fractions B(B + -> Xc\K + ) = (4.6 ± 0.4) x 10~ 4 and B(x c i -> J/ipj) = 
(34.4 ± 1.5) x 10~ 2 ; the contributions from the Xco an d Xc2 modes are neglected. 

The B + — > XcK + and B + — > J/ipK + candidates are selected keeping as many of 
the selection criteria in common as possible with the main analysis. The J/ip and Xc 
selection criteria are the same as for the prompt analysis, apart from the pseudo-decay 
time requirement. The bachelor kaon is required to have a well measured track (% 2 /ndf < 
5), a minimum impact parameter x 2 with respect to all primary vertices of greater than 
9 and a momentum greater than 5 GeV/c. The bachelor is identified as a kaon by the 
RICH detectors by requiring the difference in log-likelihoods between the kaon and pion 
hypotheses to be larger than 5. The B candidate is formed from the Xc or J/ip candidate 
and the bachelor kaon. The B vertex is required to be well measured (x 2 /ndf < 9) and 
separated from the primary vertex (flight distance x 2 > 50). The B momentum vector 
is required to point towards the primary vertex (cos# > 0.9999, where 9 is the angle 
between the B momentum and the direction between the primary and B vertices) and 
have an impact parameter x 2 smaller than 9. The combinatorial background under the Xc 
peak for the B + — > XcK + candidates is reduced by requiring the mass difference AM Xc = 
M(/i + /i~7) - M(fi + /i~) < 600MeV/c 2 . A small number of B + — s- x c K + candidates 
which form a good B + — > J/ipK + candidate are removed by requiring |M(/i + fi~^K) — 
M(jj, + trK)\ > 200MeV/c 2 . 

The AM B + = Mifi+fi-'yK) - M(/i + /i~7) mass distribution for the B + XcK + 
candidates is shown in Fig. [3|a); AM B + is computed to improve the resolution and hence 
the signal-to-background ratio. The B + — > XcK + yield, 142 ±15 candidates, is determined 
from a fit that uses a Gaussian function to describe the signal peak and a threshold 
function, 



where x = AM B + / itlq and mo, a, b and c are free parameters, to model the background. 
The reconstructed B + mass distribution for the B + — > J/ipK + candidates is shown in 



N B+ ^ XcK+ B(B + ^J/^K + ) 



x R e 



(4) 



N B+ ^ K+ B(B+ X cK+) ■ B( X c -> J/07) 




(5) 
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Fig. |3]^b). The B + — > J/vpK + yield, 8440 ± 96 candidates, is determined from a fit that 
uses a Crystal Ball function [T3] to describe the signal peak and an exponential to model 
the background. 

The photon efficiency from the observation of B + — > XcK + and B + — > J/ip K + decays 
is measured to be e 7 = (11.3 ± 1.2 ± 1.2)% where the first error is statistical and is 
dominated by the observed yield of B + — > XcK + candidates, and the second error is 
systematic and is given by the uncertainty on the branching fraction B(B + — > Xc\K + ). 
The photon efficiency measured in data can be compared to the photon efficiency, (11.7 ± 
0.3)%, obtained using the same procedure on simulated events. The measurements are 
in good agreement and the uncertainty on the difference between data and simulation 
is propagated as a ±14% relative systematic uncertainty on the photon efficiency in the 
measurement of cr(Xc — >• J/07) / cr ( t VV')- 




Figure 3: (a) Reconstructed AM B + = M (fi + \i~ 7 K) — M(/i + /i~7) mass distribution for B + — > 
XcK + candidates and (b) the reconstructed B + mass distribution for B + — > J/ipK + candidates. 
The LHCb data are shown as solid black points, the full fit functions with a solid blue (up- 
per) curve, the contribution from signal candidates with a dashed red (lower curve) and the 
background with a dashed purple curve. 



5. Polarisation 

The simulation used to calculate the efficiencies and, hence, extract the result of Eq. ([T]) 
assumes that the J/0 and Xc are unpolarised. The effect of polarised states is studied by 
reweighting the simulated events according to different polarisation scenarios; the results 
are shown in Table [Tj It is also noted that, since the ip(2S) decays predominantly to 
J/ip 7T7T, with the 7T7T in an S wave state [2], and the ip{2S) polarisation should not differ 
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significantly from the polarisation of directly produced J/0 mesons, the effect of the 
polarisation can be considered independent of the '0(25') — > J/ipX contribution [15] . 

The J/if} and Xc — > J/07 angular distributions are calculated in the helicity frame 
assuming azimuthal symmetry. This choice of reference frame provides an estimate of the 
effect of polarisation on the results, pending the direct measurements of the J/0 and Xc 
polarisations. The J/0 system is described by the angle 9j/^ , which is the angle between 
the directions of the /i + in the J/0 rest frame and the J/0 in the laboratory frame. The 
distribution depends on the parameter which describes the J/0> polarisation; 
Aj/t/, = +1,-1,0 corresponds to pure transverse, pure longitudinal and no polarisation 
respectively. The Xc — > J/07 system is described by three angles: Xc and 0, where 
is the angle between the directions of the /i + in the J/0 rest frame and the J/0 
in the Xc rest frame, # Xc is the angle between the directions of the J/0 in the Xc rest 
frame and the Xc in the laboratory frame, and is the angle between the J/0 decay 
plane in the Xc rest frame and the plane formed by the \ c direction in the laboratory 
frame and the direction of the J/0 in the Xc rest frame. The general expressions for the 
angular distributions are independent of the choice of polarisation axis (here chosen as 
the direction of the Xc in the laboratory frame) and are detailed in Ref. jl]. The angular 
distributions of the Xc states depend on m XcJ which is the azimuthal angular momentum 
quantum number of the Xcj state. 

For each simulated event in the unpolarised sample, a weight is calculated from the 
distributions of Q'm, Xc and in the various polarisation hypotheses compared to the 
unpolarised distributions. The weights shown in Table [T] are then the average of these 
per-event weights in the simulated sample. For a given (|m Xcl |, |-m Xc2 1 , Xm) polarisa- 
tion combination, the central value of the determined cross-section ratio in each pi^ bin 
should be multiplied by the number in the table. The maximum effect from the possi- 
ble polarisation of the J/0 , Xci and Xc2 mesons is given separately from the systematic 
uncertainties in Table [3] and Fig. |4j 

6. Systematic uncertainties 

The systematic uncertainties detailed below are measured by repeatedly sampling 
from the distribution of the parameter under consideration. For each sampled value, the 
cross-section ratio is calculated and the 68.3% probability interval is determined from the 
resulting distribution. 

The statistical errors from the finite number of simulated events used for the calcula- 
tion of the efficiencies are included as a systematic uncertainty in the final results. The 
uncertainty is determined by sampling the efficiencies used in Eq. [T] according to their er- 
rors. The relative systematic uncertainty due to the limited size of the simulation sample 
is found to be in the range (0.3 — 3.2)% and is given for each bin in Table [2 

The efficiency extracted from the simulation sample for reconstructing and selecting 
a photon in Xc — > J/07 decays has been validated using B + — > XcK + and B + — > J/ijjK + 
decays observed in the data, as described in Sect. |4j The relative uncertainty between 
the photon efficiencies measured in the data and simulation, ±14%, arises from the finite 
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Table 1: Polarisation weights in p T ' bins for different combinations of the J/ip, Xci aud Xc2 
polarisations. is the J/ip polarisation parameter; Aj/y, = +1,-1,0 corresponds to fully 

transverse, fully longitudinal and no polarisation respectively. m XcJ is the azimuthal angular 
momentum quantum number corresponding to total angular momentum J; Unpol means the \c 
is unpolarised. 



^ 1 Xcl 1' 1 Xe2 1' J/w ' 
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Table 2: Summary of the systematic uncertainties on a{\ c — > J/ipj) I v{J/'4>) i n eacri p4 bin. 



7 /ih 

( GeV/c) 


2-3 


3-4 


4-5 


5-6 


6-7 


7-8 


Size of simulation sample 


+0.0006 
-0.0005 


+0.0006 
-0.0005 


+0.0007 
-0.0006 


+0.0009 
-0.0009 


+0.001 
-0.001 


+0.002 
-0.002 


Photon efficiency 


+0.011 
-0.010 


+0.013 
-0.011 


+0.013 

-0.012 


+0.016 
-0.013 


+0.016 

-0.013 


+0.017 
-0.015 


Non-prompt J/0 fraction 


+0.002 
-0.005 


+0.003 
-0.005 


+0.003 
-0.006 


+0.004 
-0.008 


+0.005 
-0.010 


+0.006 
-0.011 


Fit model 


+0.003 
-0.003 


+0.003 

-0.003 


+0.002 

-0.004 


+0.003 
-0.005 


+0.002 
-0.005 


+0.003 
-0.006 


Simulation calibration 


+0.010 
-0.000 


+0.010 

-0.000 


+0.012 

-0.000 


+0.012 
-0.000 


+0.015 

-0.000 


+0.014 
-0.000 


( GcV/c) 


8-9 
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13-15 


Size of simulation sample 


+0.002 
-0.002 


+0.003 
-0.003 


+0.004 
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+0.006 
-0.006 


+0.008 
-0.008 


+0.008 
-0.008 


Photon efficiency 
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-0.016 


+0.020 
-0.016 


+0.019 

-0.016 


+0.019 
-0.018 


+0.021 
-0.020 


+0.023 
-0.019 


Non-prompt J/0 fraction 


+0.009 
-0.011 


+0.012 

-0.013 


+0.011 

-0.017 


+0.019 
-0.019 


+0.022 
-0.018 


+0.018 
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Fit model 


+0.002 
-0.005 


+0.002 
-0.003 


+0.006 
-0.002 


+0.001 
-0.006 


+0.003 

-0.008 
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Simulation calibration 
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-0.000 
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size of the observed B + — > XcK + yield and the uncertainty on the known B + — > XciK + 
branching fraction, and is taken to be the systematic error assigned to the photon efficiency 
in the measurement of <j(xc J 1^1) I <^(J/ip). The relative systematic uncertainty on the 
cross-section ratio used in Eq. [Tjis determined by sampling the photon efficiency according 
to its systematic error. It is found to be in the range (6.4 — 8.7)% and is given for each 
p^ bin in Table 2 

The J/ip yield used in Eq. [I] is corrected for the fraction of non-prompt J/ip , taken 
from Ref. pj. For those pi^ and rapidity bins used in this analysis and not covered by 
Ref. [2] (13 < p^ < 14 GeV/c and 3.5 < y J ^ < 4.5; 11 < < 13 GeV/c and 4 < y J ^ < 4.5; 
and 14 <p^ < 15GeV/c), a linear extrapolation is performed, allowing for asymmetric 
errors. The systematic uncertainty on the cross-section ratio is determined by sampling 
the non-prompt J/ip fraction according to a bifurcated Gaussian function. The relative 
systematic uncertainty from the non-prompt J/ip fraction is found to be in the range 
(1.3 — 10.7)% and is given for each p^ bin in Table 2J 

The method used to determine the systematic uncertainty due to the fit procedure in 
the extraction of the Xc yields is discussed in detail in Ref. [3]. The uncertainty includes 
contributions from uncertainties on the fixed parameters, the fit range and the shape of 
the overall fit function. The overall relative systematic uncertainty from the fit is found 
to be in the range (0.4 — 3.2)% and is given for each bin of p^ in Table 2} 

The systematic uncertainty related to the calibration of the simulation sample is eval- 
uated by performing the full analysis using simulated events and comparing to the ex- 
pected cross-section ratio from simulated signal events. The results give an underestimate 
of 10.9% in the measurement of the a(Xc — > Jfyl) / ^(J/ip) cross-section ratio. This de- 
viation is caused by non-Gaussian signal shapes in the simulation which arise from an 
untuned calorimeter calibration. These are not seen in the data, which is well described 
by Gaussian signal shapes. This deviation is included as a systematic error by sampling 
from the negative half of a Gaussian with zero mean and a width of 10.9%. The rela- 
tive uncertainty on the cross-section ratio is found to be in the range (6.3 — 8.2)% and 
is given for each bin of in Table 2 A second check of the procedure was performed 
using simulated events generated according to the distributions observed in the data, i.e. 
three overlapping Gaussians and a background shape similar to that in Fig. [TJ In this 
case no evidence for a deviation was observed. Other systematic uncertainties due to the 
modelling of the detector in the simulation are negligible. 

In summary, the overall systematic uncertainty is evaluated by simultaneously sam- 
pling the deviation of the cross-section ratio from the central value, using the distributions 
of the cross-section ratios described above. The systematic uncertainty is then determined 
from the resulting distribution as described earlier in this section. The separate system- 
atic uncertainties are shown in bins of p!^ in Table [2] and the combined uncertainties are 
shown in Table [U 
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Table 3: Ratio a(Xc ~ > Jfyl) I a (J/' l P) i n bins of pJ* i n the range 2 <p T ' < 15 GeV/c and in 
the rapidity range 2.0 <y J ^ <4.5. The first error is statistical and the second is systematic 
(apart from the polarisation). Also given is the maximum effect of the unknown polarisations 
on the results as described in Sect. HJ 



pf^GeV/c) 




Polarisation effects 


2-3 


n 1 /in+ - 005 +0-015 

U.±*±U_ .005 -0.011 


+0.025 
-0.014 


3-4 


n 1 «n+ - 003 +0-017 

U.1UU_ 004 -0.012 


+0.028 
-0.015 


4-5 


D 1 «Q+0.003 +0.019 
U.1UO_ 003 -0.012 


+0.035 
-0.018 


5-6 


D 1 8Q+ - 004 +0.021 
U.lO»_ 004 -0.015 


+0.048 
-0.025 


6-7 


D 1 8Q+ - 005 +0.022 

u.ioy_ 004 -0.016 


+0.054 
-0.028 


7-8 


n 9-1 -i +0.005 +0.024 
U.Z11_ 005 _ .017 


+0.064 
-0.033 


8-9 


D 91 O+0.007 +0.026 
U.ZlO_ 007 _ .019 


+0.068 
-0.034 


9-10 


n 990+O.OO9 +0.030 
U.ZZO_ 009 _ .019 


+0.070 
-0.034 


10 - 11 


D O9«+0.011 +0.030 
U.ZZU_ Q11 _ .022 


+0.073 
-0.036 


11 - 12 


D 9QQ+0.013 +0.034 
U.ZOO_ Q 013 _ .026 


+0.070 
-0.036 


12 - 13 


D 9^9+0.018 +0.037 
U.ZJZ_ Q17 _ .029 


+0.071 
-0.035 


13-15 


D 9«8+ - 018 +0-038 
U.ZUO_ .oi7 -0.025 


+0.080 
-0.037 



7. Results and conclusions 

The cross-section ratio, a(x c — > Jfyl) I ^(J/i^), measured in bins of p^ is given 
in Table [3] and shown in Fig. |4j The measurements are consistent with, but suggest a 
different trend to previous results from CDF using pp collisions at -y/s = 1.8 TeV [5] as 
shown in Fig. 4^a), and from HERA- 1? in pA collisions at y/s = 41.6 GeV, with p^ below 
roughly 5 GeV/c, which gave a(xc ->■ J/^i) / cr(J/ip) = 0.188 ± 0.013±g;gH @]. 

Theory predictions, calculated in the LHCb rapidity range 2.0 <y J ^ <4.5, from the 
ChiGen Monte Carlo generator [16] and from the NLO NRQCD calculations [17] are 
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Figure 4: Ratio cr(Xc — > Jfyl) / (r(J/ij>) in bins of p T in the range 2<p T / < 15 GeV/c. The 
LHCb results, in the rapidity range 2.0 < y J ^ < 4.5 and assuming the production of unpolarised 
J/tj) and Xc mesons, are shown with solid black circles and the internal error bars correspond to 
the statistical error; the external error bars include the contribution from the systematic uncer- 
tainties (apart from the polarisation). The lines surrounding the data points show the maximum 
effect of the unknown J/t/j and Xc polarisations on the result. The upper and lower limits cor- 
respond to the spin states as described in the text. The CDF data points, at y/s = 1.8 TeV in 
pp collisions and in the J/ip pseudo-rapidity range \rj J '^ \ < 1.0, are shown in (a) with open blue 
circles [5]. The two hatched bands in (b) correspond to the ChiGen Monte Carlo generator 
prediction [E] and NLO NRQCD [H]. 
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shown as hatched bands in Fig. Nib). The ChiGen Monte Carlo event generator is an 
implementation of the leading-order colour-singlet model described in Ref. [TS] . However, 
since the colour-singlet model implemented in ChiGen does not reliably predict the 
prompt J/0 cross-section, the cr(xc — » J/07) / <r(J/0) prediction uses the J/0 cross- 
section measurement from Ref. [2] as the denominator in the cross- section ratio. 

Figure [4] also shows the maximum effect of the unknown J/0 and Xc polarisations on 
the result, shown as lines surrounding the data points. In the first pi^ bin, the upper 
limit corresponds to a spin state combination (|m Xcl |, |m Xc2 |, equal to (1, 2, —1) and 
the lower limit to (0, 1, 1). For all subsequent bins, the upper and lower limits correspond 
to the spin state combinations (0, 2, —1) and (1, 0, 1) respectively. 

In summary, the ratio of the u{Xc Jfyl) / o~(J/ip) prompt production cross-sections 
is measured using 36 pb -1 of data collected by LHCb during 2010 at a centre-of-mass 
energy = 7 TeV. The results provide a significant statistical improvement compared 
to previous measurements jU The results are in agreement with the NLO NRQCD 
model [T7j over the full range of p^ . However, there is a significant discrepancy com- 
pared to the leading-order colour-singlet model described by the ChiGen Monte Carlo 
generator [16]. At high p^ , NLO corrections fall less slowly with pi^ and become im- 
portant, it is therefore not unexpected that the model lies below the data. At low p^ , 
the data appear to put a severe strain on the colour-singlet model. 
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